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Abstract 

Two different approximation schemes for the self-consistent solution of the 
relativistic Brueckner-Hartree-Fock equation for finite nuclei are discussed using 
realistic One-Boson-Exchange potentials. In a first scheme, the effects of corre- 
lations are deduced from a study of nuclear matter and parameterized in terms 
of an effective a, u and tt exchange. Employing this effective interaction rela- 
tivistic Hartree-Fock equations are solved for finite nuclei 16 0, 40 Ca and 48 Ca. 
In the second approach the effect of correlations are treated in the Brueckner- 
Hartree-Fock approximation directly for the finite nuclei, but the modifications 
of the Dirac spinors in the medium are derived from nuclear matter assum- 
ing a local-density approximation. Both approaches yield rather similar results 
for binding energies and radii in fair agreement with experimental data. The 
importance of the density dependent correlation effects is demonstrated and 
different ingredients to the spin-orbit splitting in the shell-model of the nucleus 
are discussed. 
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1 Introduction 

The various attempts to derive the bulk properties of nuclear systems from realistic 
nucleon-nucleon (NN) interactions are confronted with two major obstacles. The first 
one is the necessity to consider the effects of NN correlations which are due to the 
strong short-range and tensor components in a realistic NN interaction. The second 
one is of a relativistic nature: the strong scalar-meson (a) exchange part required 
in realistic meson-exchange potentials |l], 0, gives rise to a significant modification 
of the Dirac structure of nucleons in the nuclear medium |J. Therefore relativistic 
features should be included in the many-body theory of nuclear systems, in order to 
account for this effect. 

The importance of the NN correlations is made obvious by the fact that no binding 
energy of nuclear systems is obtained if these correlations are ignored: A Hartree-Fock 
(HF) calculation employing e.g. a realistic One-Boson- Exchange (OBE) potential ||, 
which fits NN scattering data, yields unbound nuclei. Various methods have been 
developed to include the effects of two-nucleon correlations. One possibilitiy is the so- 
called Brueckner-Hartree-Fock (BHF) approximation. In this approach one considers 
a Slater-determinant, which should be an appropriate model wave-function for the 
nuclear system to be investigated. Solving the Bethe-Goldstone equation yields an 
effective interaction, the G-matrix, which depends on the bare NN interaction and 
the model wave-function considered. The self-consistency condition of BHF now 
requires that the model-wave function, which is needed to set up the Bethe-Goldstone 
equation, is made identical to the solution of the HF equations using the G-matrix 
SIS db kind of effective interaction. 

This self-consistency problem is simplified for nuclear matter since the transla- 
tional symmetry of this infinite system requires plane waves for the single-particle 
wave functions to built up the Slater-determinant. For many years, however, the 
BHF self-consistency problem has also been solved for finite nuclei ||. |5|] . 

The inclusion of NN correlations led to a substantial improvement in the micro- 
scopic description of bulk properties of nuclei. For both nuclear matter as well as 
finite nuclei BHF calculations employing various realistic NN interactions gave re- 
sults, which were located on the so-called Coester band || [7|]. This means that either 
the calculated binding energy turned out to be too small or the calculated radii were 
too small (which corresponds to a saturation density of nuclear matter too large) as 
compared to the experimental data. Attempts have been made to improve the many- 
body approach such that results "off the Coester band", closer to the experimental 
data were obtained. Such attempts have been made within the hole-line expansion 
of the Brueckner theory |5| or using different schemes]^, [H], |TI[] . Even today it is 
not really clear if it is possible to derive the bulk properties of nuclear systems from 



realistic NN forces within a non-relativistic many-body theory [12]. 

Motivated by the success of the phenomenological o - uj model of Walecka and 
Serot H, attempts have been made to incorporate the relativistic features of this ap- 
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proach also in nuclear structure calculations which are based upon realistic NN forces. 
In this approach, one accounts for the fact that the relativistic nucleon self-energy 
in a nuclear medium is given essentially by a large attractive component, originating 
mainly from the exchange of the scalar a meson and therefore transforming like a 
scalar under a Lorentz transformation, and a repulsive component, which transforms 
like the zero component of a Lorentz vector and is mainly due to the exchange of the 
u meson. The single-particle motion is described by a Dirac equation which includes 
this self-energy. The strong components of the self-energy yield solutions of the Dirac 
equation, which are quite different from the Dirac spinors describing the nucleons in 
the vacuum. 

This change of the Dirac spinors in the medium gives rise to a self-consistency 
problem beyond the one already discussed above: Already in evaluating the matrix 
elements of the bare NN force V one should know the structure of the Dirac spinors, 
resulting from the solution of the Dirac equation. Again, this self-consistency problem 
is highly simplified in nuclear matter. In this case the medium-dependence of the 
Dirac spinors is characterized by an effective mass, which represents the ratio of the 
small to the large component of the spinor. 

Such Dirac BHF (DBHF) calculations have been performed for nuclear matter by, 
e. g., Shakin and collaborators ||13]| , Brockmann and Machleidt |TJ|, and ter Haar and 



Malfliet [15]. The basic aspects of this approach have been thoroughly investigated 



by Horowitz and Serot [0. Due to the scalar field, the nucleon mass is reduced 
enhancing the ratio between small and large components of the Dirac spinors. This 
change in the Dirac spinors yields a reduction of the scalar density, which implies 
that the attraction due to the exchange of the a meson in OBE potentials is reduced. 
At small densities of nuclear matter this loss of attraction is counterbalanced by a 
reduction of the kinetic energy, which is also caused by the medium dependence of 
the Dirac spinors. At larger densities the loss of attraction in the NN interaction 
overwhelms the loss of repulsion in the kinetic energy and for those densities the en- 
ergy calculated in the DBHF approximation is less attractive than the corresponding 
energy calculated in the BHF approximation ignoring these relativistic effects. 

Consequently the saturation points calculated for nuclear matter in DBHF approx- 
imation are shifted to smaller densities as compared to the BHF result. Brockmann 
and Machleidt succeeded in constructing a realistic OBE potential which fits NN scat- 
tering data and also yields DBHF results for nuclear matter in satisfying agreement 
with the empirical data ||14|| . The same feature is also observed for the potential " A" , 
defined in table A. 2 of ref.@, which we will consider also in our present investigation. 

This success of the DBHF approximation in nuclear matter gives rise to the hope 
that the same DBHF approximation may also be successful to reproduce the binding 
energies and radii of finite nuclei. From the discussion above, it is obvious, however, 
that a complete self-consistent calculation for finite nuclei is rather involved. There- 
fore we are going to investigate two approximations, in which either the effects of 
correlations or the relativistic effects are taken from studies of nuclear matter, while 
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the respective other components of the calculation are treated in a self-consistent way 
directly for the finite nuclei. 

In the first approximation, we determine an effective meson theory (er, uo and tt 
mesons), which yields in a Hartree-Fock approximation for nuclear matter at a given 
density p the same observables for the self-energy of the nucleons and the binding 
energy as a DBHF calculation of nuclear matter. This leads to a set of coupling 
constants, which depend on the nuclear density. This density dependence reflects the 
density dependence of the correlations described by the G-matrix of DBHF. Keeping 
track of the density dependence of these coupling constants one can perform a rela- 
tivistic HF calculation using techniques as described e.g. by Bouyssy et.al. |T7[. A 



calculation along this line has been performed by Brockmann and Toki (JT8J restricted 
to a Hartree description and first results of a HF calculation have been presented in 



l!J . Both of these earlier investigations allowed for an exchange of effective a and u 
mesons, only. 

In the second approximation one is expanding the single-particle wave functions of 
a self-consistent BHF calculation for finite nuclei in a basis of plane waves. The Dirac 
structure of these plane waves is taken from the corresponding state in nuclear matter 
of a density, which is equal to the average density for the single-particle orbit in the 
finite nucleus under consideration. In this way one deduces the Dirac effects from 
nuclear matter, but otherwise performs a complete self-consistent BHF calculation. 
This approach has also been used in ref . p(] . 

It turns out that these two very different approximations yield rather similar 
results. Therefore one can assume that a result of complete DBHF calculation should 
also be close. It is the main aim of the present investigation to explore the differences 
between the two approaches in a systematic way. For that purpose we are studying 
results on various nuclei ( 16 0, 40 Ca, 48 Ca) using different OBE interactions (OBE 
potentials "A" and "C" defined in table A. 2 of @). After this introduction we define 
the details of the two approximations towards a self-consistent DBHF calculation for 
finite nuclei in sections 2 and 3, respectively. The results are presented and discussed 
in section 4 and the main conclusions are summarized in section 5. 



2 Effective Meson Exchange Approach 

Our starting point for the description of the nuclear many - body - problem is the 
effective lagrangian density for the interacting nucleons and the a, u and 7r meson 

C = Co + &, (1) 
consisting of the free lagrangian density 

C = %3"-M)* + i(^<D/$,-m 2 X) 

+\ml^t ~ jv*" + \{%**&*« - m l€) (2) 
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with 

= d^ WtV - d v ^ Utll , (3) 
and the interaction lagrangian density 

& = -G^m^-G^^- A^ T57 (^$ ff )*. ( 4 ) 

The nucleon field and rest mass is denoted by \1/ and M, whereas the meson 
fields, rest masses and effective nucleon- meson coupling constants are denoted by $j, 
mi and Gi or with z = {a, u, tt} for the scalar, vector and pseudoscalar meson 
respectively. Note that for the pion, we use the pseudovector coupling and suppress 
the notations for the isospin degrees of freedom. Moreover we already mention here, 
that we will subtract the zero range part in the one-pion-exchange contributions to 
the NN interaction. This is done to account for the effects of short-range correlations 
between the interacting nucleons. 



2.1 Nuclear Matter 



Following standard techniques [|3], [IT]], the Hartree-Fock approximation for this la- 
grangian leads to a Dirac equation for nucleons with four-momentum p = (p , p) in 
nuclear matter 

[ 7 • p + M + £(p)]*(p, s) = 7 o£(pMp, s) , (5) 

for the nucleon spinors ^(p, s), containing the self energy S(p). Because of the 
isotropy of nuclear matter the spinors ^(p, s) are known to be plane waves and in the 
rest frame of nuclear matter the self energy E(p) for on-shell nucleons (p° = E(p)) de- 
pends only on the absolute value of the three-momentum p. This nucleon self energy 
can be split into different parts with a well-defined behavior under Lorentz transfor- 
mations. Because of parity conservation, time reversal invariance and hermiticity the 
most general form of S(p) is restricted to 

£(p) = £*(p)-7 £°(p)+7-p£>) (6) 

with S s (p), and ^ v (p) transforming like Lorentz scalars. Therefore the Dirac 

equation can be rewritten as 

[ 7 • p* + M*(p)]*(p, s) = 7 o£(p)**(p, s) , (7) 
introducing the definitions 

p* =p(l + £*(p)) , 
M*(p) =M + S s (p), 

E*(p) = E(p) + E°(p) . (8) 
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The formal similarity with a free Dirac equation allows immediately to determine the 
nucleon spinor in nuclear matter to 



*(p,s) 



£*(p) + M*(p) 
. 2E*(p) 



1/2 



1 



E*(p)+M*(p) 



(9) 



now with a modified ratio of the spinors upper and lower component if compared 
with the vacuum solution. The spinors are normalized (noncovariant) to 



#(p,s)%(p,s) = 1 , *(p,s)*(p,s) 



M*(p) 



E*(p) 



and the on-shell condition in nuclear matter now reads like 

\2 T\/T*(„\2 i „*2 



E*(pY = m*(pY + P * 



(10) 



(11) 



On the level of the Hartree-Fock approximation, the mesons used in our lagrange 
density gives rise to the following contributions to the self energy: 



s s (p) 



+ 



l l r k F 



(Air) 2 p Jo 



qdq 



M*(q) 



GlQ a (p,q)-4GlQ^p,q) 



E°(p) = -p) 



- 3 ( — J mlQ^p.q) 
G u 



(12) 



1 1 



(47r) 2 p Jo 



qdq 



GlQ tr (p,q) + 2G 2 w Q w (p,q) 



-3 



E"(P) 



( u_ 

1 



mlQ v (p,q) 



(13) 



(Anp)'- 



E*(q) 



2G 2 a r a (p,q) + 4Glr u} (p,q) 



-6 



/ 

JO 



((P 2 + g 2 )I\r(p, g) - pq®v(p, q 



(14) 



We omit in our notation the obvious dependence of the self energies on the Fermi 
momentum A; jr. The first term in X s (p) and £°(p) corresponds to the Hartree con- 
tribution using 



p(k F ) 



3tt 2 



k 3 F and p s (k F ) = — / g 2 dg 



7T 2 Jo 



(15) 
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for the baryon and scalar density, respectively. The remaining expressions are due to 
the Fock (exchange) contributions where we have used the abbreviations 

A(p,q) = p 2 + q 2 + m 2 -(E(p)-E(q)) 2 (16) 
rw n , ( Ai(jp,q) + 2pq\ 

B ' M = '" Um-W (17) 

r ((p , 5 ) = ^^'' -1, (is) 

again z = {a, u;,7r}. Two important things have to be noted. First, as already 
mentioned above, we have subtracted zero-range contributions from pion exchange to 
the self energy £(p). Secondly, we don't want to consider retardation effects in the 
meson propagators. Although retardation causes no problems in nuclear matter, the 
neglection leads to significant simplifications in finite nuclei. 

In the next step we determine effective coupling constants G a and for the scalar 
and vector meson by requesting that the HF expressions for the scalar self energy 
S s (p) calculated at the Fermi surface (p = kf) and the binding energy per nucleon 
reproduce the corresponding results of a Dirac-Brueckner-HF (DBHF) calculation 
[0, [TJJ using realistic NN-forces, namely versions A and C of the Bonn potential 
0]. For the Pion we fix the coupling constant to the free value f 2 /4:ir = 0.08 and the 
masses of the mesons are chosen to be identical to those of the OBE potential (m 7r =138 
MeV, m cr =550 MeV and m w =783 MeV). In this way we obtain for each baryon 
density p two effective coupling constants G a (p) and G u {p). The density dependence 
of these coupling constants reflects the density dependence of the correlations taken 
into account in the DBHF approximation. 



2.2 Finite Nuclei 

Once the density dependent coupling constants are determined, we assume for the 
study of finite nuclei, that we can account for the density dependent correlation effects 
in a relativistic HF calculation by employing the coupling constants calculated at the 
local density G a (p(r)) and G u {p{r)), where the density profile p(r) is determined 
from the result of the relativistic HF calculation in a self-consistent manner. For that 
purpose we write the nucleon spinor for the finite system in coordinate space 

< r\a >= * Q (f) = (_J^^^y Ka , m Mx Uqa) (19) 
_ ( g a (r) y Ka ,m a (ty \ , s 

- {if a (r)y- Ka , m M XUqa) - (20) 



Again, the spinors are normalized to 

d 3 r$?Hr)iS?(r) = I r 2 dr 



o 



9l(r) + f 2 a (r) =1- (21) 
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All quantum numbers are summarized by the index a = {a, m a } with a = {n a , n a , q a }- 
n a characterizes the radial quantum numbers, whereas n a = (2j a +l)(l a —j a ) describes 
the angular momenta. Obviously the upper and lower spinor-component for the same 
total angular momentum j a have different orbital quantum numbers l a . We introduce 
the corresponding l' a to the same j a 

]' _ J ~l~ 1 ^ or la — ja ~ 2 (00\ 

La ~\ / a -lfor/ a =j a + ± [ZZ) 
and a' = {n a , n' a , q a }, n' a = (2j a + l)(l' a - j a ). y Ka , ma (ty is constructed as usual 

y Ka ,m a (ty = J2 ( l a m l l / 2m s\ja m a)Yl a , mi (tyXl,m s ■ ( 2 3) 
mi i m s 

For the isospinor X\,(q a ) we use 

J +1 for protons , , 

^ a 1 — 1 for neutrons 

The Dirac equation is solved by expanding the radial functions g a (r) and f a (r) in a 
discrete basis of spherical Bessel functions. The wave numbers for this basis are chosen 
such that this discrete basis is a complete orthonormal basis in a sphere of radius D. 
This radius is chosen to be large enough that the results for the bound single-particle 
states are independent on D. With this expansion the Dirac equation (^) is rewritten 
in form of an eigenvalue problem and the eigenvalues (E a ) and eigenvectors (the 
expansion coefficients for g a and f a ) are determined by a simple matrix diagonalisation 

In the following we give the expressions for the matrixelements of the self-energy 
in this Dirac matrix, calculated in the Hartree-Fock approximation. In this work we 
consider nuclei with a closed proton and neutron shell only. Therefore, the isovec- 
tor pseudoscalar meson yields no contributions in the Hartree approximation. The 
Hartree matrix elements for the isoscalar scalar and vector part of the interaction and 
the Coulomb force are given as 

/•OO 

<a|£f|/3> = S Ka , Kb 5 ma>mb 8 qajqb r 2 drG a {r)[g a {r)g b {r) - f a (r)f b (r)] 

Jo 

x J r n drV a {r')p s {r')lQ{m a r < )%{m a r > )\ (25) 

/•OO 

< «|£ff|/3 > = S KatKb 5 ma>mb 5 qajqb / r 2 drG w (r) [g a (r)g b (r) + f a (r)f b (r)\ 



o 



x\m u r /2 dr'G (J (r')p(r')/o(m a; r < )Aro(m w r > )^ (26) 



/■OO 

<a|£f|/3> = 5 Ka>Kb S ma>mb r 2 dr[g a (r)g b (r) + f a (r)f b (r)} 



o 



x L.fy,fr/ftWlfl + fc>fl+») (27) 
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with the definitions 



p(r) 
Ps{r) 



a 

^Eia [gl{r)-Urf 



{21 



for the baryon and scalar density in finite nuclei, where j = y/2j + 1 and index a 
running over all occupied orbits. If the baryon density of the protons p p (r) is needed 
as in the Coulomb self energy , index a runs only over occupied proton orbits. 
Note the use of the coupling constants depending on the local position, which is just 
an abbreviation for Gi{r) = Gi(p(r)). In eqs.(p5|) - r< ancFr> are denoting the 
smaller or bigger value of r and r' . The functions Il{x) and Kl(x) arise from the 
multipole expansion of the meson propagator in coordinate space and are defined 
using the modified spherical Bessel functions I and i<T[25|: 



h(x) = h+l/ ^ X) K L (x) = KL+1 £ {X) 
\/x Jx 



(29) 



The Fock contributions originating from the a exchange read 



-m. 



< &\K\f3 >= $Ka,K b $m a ,m b 5 qa , qb K 

J a c 



Qa,qc 



x / r 2 dr\ G a (r) [g a (r)g c (r) - f a {r)f c {r)} 



x / r' 2 dr'G a (r') [g c {r')g h {r') - f c {r')f b {r')] 
Jo 

x < a ll^l|c > 2 /L(m CT r < )K" L (m CT r > ) 

L 



(30) 



with index c running over all occupied states. Here and in the following expressions 
for the other mesons the < a||Y^||6 > represent the reduced matrix elements of the 
spherical harmonics and 



<a||T J (L)||6>= J^(-Y<hjaJ b Lj( Q a q £ 




(31) 



2 2 



Using this notation, the matrix elements of the Fock contributions arising from the 
uj exchange can be written 



< 0^1511^1/3 > — 5 Ka jKb 5 rna:rnb 5q a: q b ^ ^q a ,qc 

J a c 

xjjf PdrG^r) [g a {r)g c {r) + /„(r)/ c (r)] 
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x / r' 2 dr'GUr')[g c (r')g b (r>) + f c (r')f b (r')] 
Jo 

< a\\Y L \\c > 2 TUm^KLim^)] 

L ' 

+ El r r2drG ^a{r)fc{r) < c\\Tj(L)\\a > 
-f a (r)g c (r)<c\\Tj(L)\\a' >] 

roo 

x / r' 2 dr'GUr')[g c (r')f b (r') < c\\Tj(L)\\a' > 
Jo 

-fc(r')g b (r') < c'\\Tj(L)\\a >]7 L (m a; r < )X L (m w r > ) \ (32) 
The isovector pseudoscalar meson using the pseudovector coupling yields 



< \P >= 5 Ka:Kb 5 ma ,m b 5 qa , qb ( ^- ) ^ E(2 - £g ,,J 

\ '' 4 7r / J a c 

x r r 2 d J_|I [0 a (r)<7 c (r) + / (r)/ c (r)] [<? 6 (r)«? c (r) + / 6 (r)/ c (r)] 
+m^L- 4 | < allelic' > | 2 

L L1L2 

X [(/C a + K c + /l(-Ll)) g a (r)g c (r) ~ (« a + «c - M-^l)) /a( r )/c(^)] 
/■oo 

x / r /2 dr' [(k 6 + « c + h(L 2 ))g b (r)g c (r) - (« 6 + k c - h(L 2 ))/ b (r)/ c (r)] 
x R(L u L 2 ,r,r') 

+ W{< a\\Tj(L)\\c> g a (r)g c (r)+ < a'\\Tj(L)\\c' > / B (r)/ C (r)] 

d LJ 

x [< aHTj^Ho ^(r)^ c (r)+ < a'HTj^Hc' > / 6 (r)/ c (r)]} (33) 

with Lj = {L — 1, L + 1} and the auxiliary functions 

R(L 1 ,L 2 ,r,r') = 9(r' - r)T Ll (m ps r)K L2 (m ps r') + 9(r - r')T L2 (m ps r')K Ll (m ps r) (35) 

using the step function 9(x — y). The expressions for the self energy show that the 
Hartree contributions can be rewritten by defining a local potential, e.g. 

roo 

<a\X*\P> = S Ka:Kb S ma:mb S qatqb / r 2 dr [g a (r)g b (r) - f a (r)f b (r)] V a (r) 

J 

= ^, K ^ ma , m6 5 (?a , (?6 /rf 3 rv[;t i ( r ) 7 0v;(r)M/ /3 (r) (36) 

with r oo _ 

V a (r) = -G a (r)m a / r ,2 rfr , p s (r')G (7 (r')/o(m (7 r < )^o(m (7 r > ) , (37) 
Jo 
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whereas the Fock contributions are obviously non-local. 

Solving the Dirac equation (|5|) with the technique mentioned above in a self- 
consistent way one can finally determine the binding energy of the nucleus with A 
nucleons as 

E= l - E a + T a -AM, (38) 

a,( OCC) 

where E a denotes the single-particle energy obtained by solving the Dirac equation 
and T a is the corresponding kinetic energy. 



3 Local Density Approximation 

In contrast to the effective meson exchange approach, which, as discussed in the 
preceding section, determines the effects of correlations from studies of nuclear matter, 
we are now considering an approximation, in which the relativistic effects are deduced 
from nuclear matter. For that purpose we consider as an example the expansion of 
an harmonic oscillator (h.o.) state in terms of plane wave spinors 

^n,j,l,m(P,p) = 



e*(p,p) + m*(p,p) / i 3Wn P ) \ 

\ 9F*fn n) V,, (Q ) ® n ^ P > (39) 

\ {P,P) \ [E*(p,p)+M*(p,p)] yjl'mKMp) ) 

where $ nj ;(p) are the momentum space h.o. wavefunctions. The structure of the 
plane wave Dirac spinor, in particular the ratio of the small to large component is 
determined by the quantities p*, M* and E* as defined in eq.(||). The values actually 
used for these quantities are taken from the DBHF calculations at a density p for 
the realistic NN interaction under consideration. In this sense the Dirac structure of 
the harmonic oscillator state defined in eq.QBTf) is derived from nuclear matter of a 
given density p. For the Dirac spinors as presented in eq.(^) one can calculate the 
matrix elements of the OBE potential under consideration employing the conventional 
techniques and identify the resulting numbers with matrix elements of a two-body 
interaction V between non-relativistic harmonic oscillator states 

< a(3\V(p)\j5 > (40) 

where a ... 8 refer to the quantum numbers of the various h.o. states and the param- 
eter p is kept to memorize that the value of this matrix element depends on a density 
parameter p, which determines the Dirac structure of the spinors used to calculate 
the matrix element. This scheme can of course be generalized to single-particle wave 
functions different from h.o. functions. In a corresponding way one can also evaluate 
the matrix elements for the operator of the kinetic energy 

tap(p) = I d 3 p^ a [ 7 • P + M] * p - M5 a /3 

= j d 3 p $* a [ MM^+nm _ M ] ^ (41) 
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For the interaction defined by the matrix elements of eq. (fH)]) one may now solve the 
Bethe-Goldstone equation 

G(Z, p) = V(p) + V(p) z _q HqQ G(Z, p) (42) 



using the standard techniques of non-relativistic BHF calculations for finite nuclei p3 
The Pauli operator Q in this equation is defined in terms of h.o. states appropriate for 
the nucleus under consideration. Beside the usual dependence on the starting energy 
Z, the matrix elements of G also depend on the density parameter p characterizing 
the structure of the Dirac spinors involved. Keeping track of this additional density 
dependence one can expand the BHF single-particle states \i > and \j > in the basis 
of h.o. states \at > 

\i>=^2c ia \a> (43) 

a 

and the expansion coefficients are determined from the solution of the eigenvalue 
problem 



E 

13 



tap(Pi)+ < a 3\ G ( Z = e i + e jiPij)\f3j > 

j,( occ) 



Ci/3 = tiCia (44) 



If one ignores in this equation the medium dependence of the Dirac spinors by putting 
Pi and pij equal to zero this eq. fl44T) together with the Bethe-Goldstone eq.([42D defines 
the conventional BHF approach for realistic OBE potentials [P4}| . in addition to the 



self-consistency requirements of this BHF approach, we now want to account for the 
medium dependence of the Dirac spinors and define an average density for nucleons 
in the orbit % by 

Pi = [ d 3 r 4>*{r)p D BHF{r)(j)i{r) (45) 



with 4>i(r) the DBHF single-particle wave function and Pdbhf{t) the radial shape 
of the baryon density obtained from this calculation. This average single-particle 
density enters into the calculation of the kinetic energy and it is also used to define 
the average density for an interacting pair of nucleons by 

Pij = \fpiP~3 ■ ( 46 ) 



4 Results and Discussion 

As a first step towards the application of the effective meson exchange approach dis- 
cussed in section 2, we have to determine the coupling constants G a (p) and G UJ (p), 
depending on the baryon density p. As discussed in subsection 2.1 this is done by 
adjusting these parameters in such a way that a mean-field calculation reproduces 
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at each density the results for the scalar self-energy X s and the binding energy per 
nucleon obtained in DBHF calculations for nuclear matter p5]. For the mean-field 
calculation we consider three different approximations. In the first approach, we just 
consider the Hartree-contributions to the self-energy and total energy (see eqs.(|l"3"D 
and ([14]) ). In the second approach we consider a o-uj model in Hartree-Fock ap- 
proximation, i.e. we keep the Hartree and the Fock terms in eqs.([L3|) to (|14D which 
originate from the exchange of a o or uj meson. This approach will be called HF(cr, to) 
or HF1. In the approach HF2 or HF(a, uj,tt) we furthermore consider the effects of 
the pion-exchange, which means that we are considering the complete set outlined in 
section 2. 

Results for the effective coupling constants determined from DBHF calculations 
for OBE potentials A and C [Q, P^fl , are listed in table 1 for various densities. For all 
three approaches considered, the effective coupling constants G a and G u decrease with 
increasing density. This is also displayed in fig. III. The decrease reflects the fact that 
the terms in the G-matrix, which are of second and higher order in the interaction, 
contain contributions, which are simulated by the exchange of a scalar and a vector 
meson 



26| . Due to the Pauli operator in the Bethe-Goldstone equation ( [42| ) and due 
to the change in the energy dominator, these contributions of higher order in the bare 
interaction V are quenched with increasing density. The importance of the density- 
dependent correlation effects parameterized in terms of these coupling constants is 
reflected by the fact that the square of the coupling constants are quenched by a 
factor 2, if the nuclear density is increased from 0.2 p to 1.4 p (po denoting the 
empirical saturation density). 

Such a substantial density dependence of the effective coupling constants of course 
effects the nuclear structure calculations. This is displayed in fig. ||, where the results 
for the binding energy per nucleon and the effective mass of the nucleon determined in 
a DBHF calculation are compared to the corresponding quantities obtained in Hartree 
calculations using effective coupling constants which are determined to reproduce the 
DBHF results at a small density (dashed curve) or large density (dotted line). It is 
obvious that the large coupling constant G a determined for small densities predicts an 
effective mass at larger density, which is considerably smaller than the one obtained 
in the self-consistent DBHF calculation. This is an example to demonstrate that 
mean field calculations employing constant coupling constants tend to overestimate 
the change of the Dirac spinors in nuclear matter at high densities. 

The density dependent coupling constants also effect the predictions for the sat- 
uration properties of nuclear matter. Note that the energy versus density curve 
obtained from the Hartree calculation with coupling constant appropriate for small 
densities yields a minimum at a low density and overestimates the binding energy 
by about 8 MeV per nucleon. The Hartree calculations using coupling constants for 
large densities (dotted line in fig. Q) predict a minimum at higher densities and a 
lower energy per nucleon. Therefore the minima of the energy versus density curves 
obtained from coupling constants determined at various densities form a "band" of 
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saturation points, which is perpendicular to the "Coester band" . 

The same phenomenon can also be observed in calculations of finite nuclei. This is 
demonstrated in table where we show results of relativistic Hartree-Fock calculations 
considering coupling constants as derived from various densities in nuclear matter. 
The binding energies displayed in this table and all subsequent ones are corrected 
for spurious center of mass effects, assuming an harmonic oscillator model, and the 
radius of the charge distribution has been evaluated from the proton density, assuming 
a radius of 0.8 fm for the charge radius of the proton. 

Inspecting the results displayed in table || one observes a sensitive dependence of 
the calculated binding energy on the choice for the coupling constants. This demon- 
strates again the importance of the density dependent correlation effects contained in 
the G-matrix of nuclear matter. It is striking to see that the calculation which yields 
the largest binding energy (p = 0.2p ) also predicts the largest radius of the charge 
distribution. So we find also for finite nuclei that the ground-state properties of fi- 
nite nuclei calculated with meson exchange parameters derived from various densities 
form a "band" which is perpendicular to the normal "Coester band" |24]. 



Furthermore its worth noting that the spin-orbit splitting deduced from the dif- 
ference in the single-particle energies e p3 / 2 and e p i/ 2 is of course largest for that in- 
teraction which yields the smallest effective mass in nuclear matter (p = 0.2p ). The 
comparison displayed in table demonstrates the importance of relativistic effects for 
the spin-orbit splitting in the nuclear shell-model |27| . 

Table [| also contains a first result which is obtained when we consider an effec- 
tive meson exchange with local coupling constants, depending on the position of the 
interacting nucleons as discussed in section 2.2. As to be expected, one finds that 
the single-particle energy for the py 2 state obtained in this self-consistent calculation 
is closer to the one obtained p = 0.2p , while the single-particle energy of the deep 
lying S1/2 state is closer to the one obtained at larger densities. It should be men- 
tioned that effective coupling constants can safely be derived from nuclear matter 
calculations only for densities as low as ~ 0.2 times the saturation density po- For 
smaller densities the conventional tools to evaluate BHF energies yield unstable re- 



sults [|TTJ. Therefore we extrapolate the coupling constants to smaller densities using 
spline functions in terms of the density. 

The next question we want to investigate is the sensitivity of the effective me- 
son exchange model on the mesons taken into account. For that purpose we have 
performed Hartree calculations for finite nuclei using the local coupling constants as 
derived from the Hartree calculations of nuclear matter (column "Hart" in table p. 
In the same way we also perform HF calculations for finite nuclei, taking into ac- 
count the effects of a and to exchange using the density-dependent coupling constants 
derived from nuclear matter ("HF1") and finally consider the complete model with 
inclusion of the pion discussed in section 2 ("HF2"). 

Results obtained for these three effective meson exchange models for the nuclei 
16 0, 40 Ca, and 48 Ca are displayed in tables |3|, £| and |5], respectively. Using the OBE 
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potential A, which yields a correct description of the saturation point of nuclear 
matter, the Hartree approximation shows fair agreement with the experimental data 
for the binding energy and radius of all three nuclei considered. Both the results 
for the radius and the binding energy are slightly below the experimental values. 
Employing potential C yields larger radii but smaller binding energies. This is the 
typical feature for 2 phase-shift equivalent potentials, the results change along the 
"Coester band" . 

The inclusion of the Fock terms in HF1 reduces the calculated radii to a significant 
extent, predicting the same or a slightly smaller binding energy as compared to the 
Hartree approach. Therefore the agreement with experiment gets worse. Furthermore 
we note that the Fock terms tend to enhance the spin-orbit splitting in the single- 
particle energies, which again deteriorates the agreement with the splitting observed 
in the experimental data. 

The pion-exchange terms included in the HF2 approximation slightly improve the 
agreement between calculation and experiment. The spin-orbit splitting is reduced 
and the binding energies are larger but the results for the radii are essentially the same 
as in the HF1 approximation. Keeping in mind the sensitivity of the calculated values 
on the density dependence of the effective coupling constants displayed in table || and 
discussed above, one may conclude, however, that all three models lead to results, 
which are rather similar. 

The main purpose of this study is to compare the predictions of the effective meson 
exchange approach to the results obtained in BHF calculations, which treat the change 
of the Dirac spinors in the medium in a local density approximation (see section 3 and 
ref-POf). For that purpose the tables [| and |] show results of this approach (identified 
as DBHF) and allow a comparison with conventional BHF calculations, which ignore 
medium dependence of the Dirac spinors completely. The differences between BHF 
and DBHF results are by far not as large as those displayed in table which reflect 
the density dependence of the correlations. Therefore we conclude that the bulk 
properties of nuclei are more sensitive to the density dependence of the correlations 
than to the medium dependence of the Dirac spinors. That is why we consider the 
approach, which treats the Dirac effects in a local density approximation (DBHF), 
to give more reliable predictions for a complete Dirac - Brueckner - Hartree - Fock 
calculation than the effective meson exchange approach, which derives the correlation 
effects from nuclear matter. 

For the case of the nucleus 16 O it has already been shown in |2(J that the inclusion 
of Dirac effects in DBHF leads to larger radii and binding energies as compared to 
the predictions of conventional BHF calculations. Thereby the agreement of the 
theoretical predictions are substantially improved. This observation is supported by 
the results on 40 Ca shown in table [|. Furthermore one finds that the DBHF results 
are in fair agreement with those obtained in the relativistic HF approximation using 
effective meson exchange. This agreement supports the conclusion that both types of 
approaches are reliable approximation for a complete Dirac Brueckner calculation. 
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Finally, we want to investigate the basic assumption of the DBHF approach which 
assumes that the Dirac spinors for the single-particle states in finite nuclei can be 
described in terms of plane wave spinors of nuclear matter. For that purpose we 
consider as an example the radial functions g a (r) and f a { r ) for the large and small 
component of the 0si/ 2 Dirac spinor calculated in a relativistic HF approach (HF1, 
OBEPA) for 40 Ca (see solid lines in fig. |3|). For the comparison we consider a Dirac 
spinor expanded in terms of spinors for nuclear matter assuming a harmonic oscillator 
expansion as in eq.(^). If we consider plane- wave Dirac spinors of the vacuum 
(kp=0, dotted line) the lower or small component is considerably weaker than the 
one resulting from the relativistic HF calculation. For the appropriate average density 
Pi as defined in eq.(|45"D the enhancement of the small component in the medium is 
fairly well described. This demonstrates that the 2 approaches not only lead to very 
similar results for the global observables like binding energy and radius, but also 
provide similar predictions for the components of the Dirac spinors. 



5 Conclusions 

Two different steps towards a self-consistent Dirac Brueckner calculation for finite 
nuclei are presented and discussed. In the effective meson exchange approximation one 
solves the relativistic Hartree-Fock equations directly for the finite system and deduces 
the correlation effects from nuclear matter. This is done in various models to study 
the importance of Fock exchange effects and the impact of the pion exchange. The 
density dependence of the effective coupling constants reflects the density dependence 
of the correlations encountered in the Brueckner G matrix. 

In an alternative approach (DBHF) the correlation effects are treated directly for 
the finite nuclei but the change of the Dirac spinors is determined from nuclear matter. 
It is demonstrated that the bulk properties of nuclei (binding energy and radius) are 
less affected by the change of the Dirac spinors than by the density dependence of 
the correlations. This implies that the approach which treats the correlation effects 
without approximation should provide more reliable results than the effective meson 
exchange approach. 

It turns out that both approximations yield very similar results. For the realistic 
OBE potential A defined in || binding energies per nucleon are obtained for 16 0, 40 Ca 
and 48 Ca, which are close to the experimental value (± 0.5 MeV). The predictions 
for the radii are still significantly below the experimental data (typically 0.2 fm). 
This might be improved by including correlations beyond the lowest order Brueckner 
theory. Recently it has been demonstrated that the inclusion of hole-hole scattering 
terms within a self-consistent Green function approach tends to improve BHF results 
in this direction p8| . 
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Table 1: DBHF results for nuclear matter derived from OBE potential A for the 
scalar part of nucleon self-energy and the binding energy per nucleon (E/A) for 
various Fermi momenta kf. The column G a and G u show the coupling constants 
which are needed to reproduce these results in a HF calculation. For this purpose 
three different models are considered: the Hartree approximation ignoring the all Fock 
contributions to the self-energy and binding energy (model "Hart"), the Hartree- Fock 
approximation ignoring the contribution of the pion-exchange (model "HFl") and the 
full model defined in section 2.1 (model "HF2"). All energies are listed in MeV and 
the Fermi momenta in unit fm -1 . For a comparison the lowest part of the table also 
shows results obtained for OBE potential C at one specific density. 





E/A 




Model 




G w 


0.80 


-7.27 


-134.3 


Hart 
HFl 
HF2 


12.436 
11.411 
11.227 


15.403 
12.941 
13.179 


1.00 


-10.62 


-209.8 


Hart 
HFl 
HF2 


11.177 
10.265 
10.104 


13.807 
11.668 
11.885 


1.20 


-13.44 


-288.8 


Hart 
HFl 
HF2 


10.059 
9.267 
9.118 


12.322 
10.470 
10.674 


1.40 


-15.59 


-374.9 


Hart 
HFl 
HF2 


9.224 
8.531 
8.389 


11.168 
9.539 
9.733 


1.50 


-14.88 


-416.3 


Hart 
HFl 
HF2 


8.851 
8.197 
8.048 


10.673 
9.145 
9.336 






OBE Potential C 




1.20 


-11.57 


-292.8 


Hart 
HFl 
HF2 


10.130 
9.297 
9.149 


12.534 
10.669 
10.869 
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Table 2: Results of relativistic HF calculations on 16 O, considering the exchange of 
effective a, uj and n mesons. The coupling constants are determined to reproduce 
DBHF results for nuclear matter (OBE potential A) at various densities: p = 0.2p 
(kp=0.8 fm -1 ), p = 0.5po (kp—1.1 fm -1 ) and p = 1.4p {kp—l.b fm -1 ). These results 
can be compared to those of a self-consistent calculation (last column) considering 
local coupling constants as discussed in section 2.2. Results are presented for the 
single-particle energies of proton states, the binding energy per nucleon (E/A, cor- 
rected for cm effects) and the radius of the charge distribution (Rch)- 





p = 0.2p 


P = 0.5p 


P = l-4p 


self- cons. 


e sl/2 [MeV] 


-54.4 


-48.8 


-46.6 


-47.1 


e p3 / 2 [MeV] 


-33.9 


-25.4 


-20.1 


-23.8 


e pl/2 [MeV] 


-18.6 


-13.7 


-11.3 


-17.7 


E/A [MeV] 


-15.17 


-9.29 


-5.38 


-7.73 


Rch [fm] 


2.52 


2.46 


2.37 


2.48 
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Table 3: Results of relativistic HF calculations on 16 O, considering various models for 
the effective meson exchange (Hart, HF1, HF2, see table 1) are compared to results of 
conventional BHF calculations and BHF calculations which account for Dirac effects 
in the way described in section 3 (DBHF). Further information see table 2. 





Hart. 


HF1 


HF2 


BHF 


DBHF 


Exp 








Potential A 






f i /o fMeVl 
e sl/2 [J-vievj 

e p3/2 [MeV] 

e pl/2 [MeV] 


-44 n 

tit:, w 

-21.5 
-15.8 


-44 

TT: . W 

-23.4 
-15.8 


-47.1 
-23.8 
-17.7 


-56.6 
-25.7 
-17.4 


-49 8 
-23.0 
-13.1 


-4D+8 
-18.4 
-12.1 


E/A [MeV] 
Rch [fm] 


-7.20 
2.57 


-7.23 
2.48 


-7.73 
2.48 


-5.95 
2.31 


-7.56 
2.46 


-7.98 
2.70 








Potential C 






e iV2 [MeV] 
e p3 / 2 [MeV] 
e pl/2 [MeV] 


-37.0 
-17.7 
-13.3 


-37.4 
-19.4 
-13.5 


-40.2 
-19.7 
-14.9 


-45.2 
-19.5 
-13.7 


-40.9 
-18.0 
-11.0 


-40±8 
-18.4 
-12.1 


E/A [MeV] 
Rch [fm] 


-5.60 
2.73 


-5.59 
2.62 


-6.09 
2.62 


-4.03 
2.48 


-5.30 
2.59 


-7.98 
2.70 
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Table 4: Results of relativistic HF calculations on 40 Ca. Further information see table 
3. 





Hart. 


HF1 


HF2 


BHF 


DBHF 


Exp 








Potential A 






e d5/2 [MeV] 
ei c i/9 [Me VI 
e d 3/2 [MeV] 


-19.8 
-15.4 
-14.5 


-21.0 
-13.7 
-13.2 


-20.8 
-14.1 
-14.2 


-30.2 
-24.5 
-16.5 


-21.9 
-13.8 
-10.2 


-14±2 
-10±1 
-7±1 


E/A [MeV] 

R C h M 


-8.21 
3.35 


-7.76 
3.14 


-8.09 
3.14 


-8.29 
2.64 


-8.64 
3.05 


-8.50 
3.50 








Potential C 






[MeV] 
e lsl/2 [MeV] 
e d3 /2 [MeV] 


-15.3 
-10.9 
-10.5 


-16.9 
-11.3 
-10.8 


-16.7 
-11.5 
-11.5 


-21.0 
-16.9 
-12.0 


-16.5 
-10.6 
-8.0 


-14±2 
-10±1 
-7±1 


E/A [MeV] 
[fm] 


-5.83 
3.44 


-5.80 
3.31 


-6.14 
3.32 


-5.06 
2.87 


-5.91 
3.21 


-8.55 
3.50 
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Table 5: Results of relativistic HF calculations on 48 Ca using the OBE potential A. 
Further information see table 3. 





Hart. 


HF1 


HF2 


Exp 


e d5/2 [MeV] 
ei iV2 [MeV] 
e d3/2 [MeV] 


-24.6 
-18.7 
-19.5 


-29.0 
-19.5 
-21.6 


-27.2 
-20.2 
-25.1 


-20±1 
-15.8 
-15.3 


E/A [MeV] 
Rch [fm] 


-8.35 
3.34 


-7.83 
3.15 


-7.90 
3.16 


-8.70 
3.50 
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Figure 1: Effective coupling constants for the exchange of a scalar meson (G a ) and 
vector meson (G u ) as a function of the density. The coupling constants are derived 
from DBHF calculations of nuclear matter employing the OBE potentials A (left part) 
and C (right part). Various mean field approximations are considered: the Hartree 
approximation, the Hartree-Fock approximation with only a and uo exchange (HF1) 
and the complete model of section 2, also including the pion (HF2) 
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Figure 2: Results for the binding energy per nucleon (upper part) and the effective 
mass M* of the nucleon as a function of the Fermi momentum in nuclear matter. Re- 
sults obtained from DBHF calculations (OBE potential A, solid line) are compared 
to those of Hartree calculations in which the effective coupling constants are deter- 
mined to reproduce the DBHF results at p = 0.2p (kp—0.8, long dashed line) and 
at p = l.lpo (k F =lA0, dotted line). 
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Figure 3: The radial function for the large upper component g(r) and the small 
lower component f(r) (both multiplied by r) for the Osi/2 Dirac spinor obtained in 
a relativistic HF calculation (HF1, OBEPA) for 40 Ca are compared to a harmonic 
oscillator spinor of eq.(39) assuming no medium correction (kp = 0, dotted line) or 
as predicted by nuclear matter at an average density (dashed line). 
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